Statistics and Geometrical Picture of Ring Polymer Melts and Solutions 
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We present a detailed account of a recently proposed phenomenological theory for noncatenated 
ring polymer melts (Phys. Rev. Lett. 106, 167802 (2011)). A basic assumption lies in the 
implementation of the noncatenation constraint via the effective excluded-volume effect, from which 
a geometrical picture of melts emerges. The result captures many of the salient features observed so 
far, including (i) the overall spacial size of rings, (ii) the coordinate number, i.e., the number of rings 
surrounding a given ring, (iii) the topological length scale as a function of the molecular weight and 
(iv) the effect of the chain stiffness and concentration. We also suggest a geometrical interpretation 
of the topological length scale, which may shed some light on the entanglement concept in polymeric 
systems. 

PACS numbers: 61.25.H-,36.20.-r,83.80.Sg 



I. INTRODUCTION 

The phenomenon of entanglement is ubiquitous in sys- 
tems of long flexible ropes and, once formed, it often 
has crucial consequences on the subsequent operations, 
as we know from our daily experience It also has 
a fundamental importance in molecular systems, dictat- 
ing the dynamics and rheology of dense linear polymer 
solutions of high molecular weig ht @-[3. Entanglement 
originates, in any case, from the topological constraints 
imposed by the non-crossability of the chain. In the his- 
tory of polymer science, predicting the macroscopic vis- 
coclasticity of the entangled polymer solution from the 
molecular theory has been one of the major goals. Our 
current understanding on this issue owes to the introduc- 
tion of the following two theoretical models. The tube 
model, originally proposed in the theory of rubber elas- 
ticity, assumes that the topological constraint imposed by 
the non-crossability of the chain effectively confines each 
chain inside a tube-like region, thereby reducing the orig- 
inal many-chain problem to single-chain problem in the 
"mean-field" Q . The reptation model then describes the 
overall motion of the chain under the temporal topolog- 
ical constraint along the primitive path, where the pres- 
ence of the chain end is essential for the disengaging from 
the original tube and the subsequent tube renewal @ . 

In this context, an interesting question arises in the 
system of ring polymers, also called circular or cyclic 
polymers 0, H[. Here, by definition, there are no chain 
ends, and the topology of the system is frozen at the 
stage of the sample synthesis, the consequence of which 
is pronounced even for static properties. To emphasize 
the point, let us recall the basic static properties of linear 
polymers [H 0, Q . In dilute solutions, their conformation 
can be characterized as a self-avoiding walk with the over- 
all chain dimension Rq ~ aN V!i , where a is the segment 



size and N is the number of segments per chain. Un- 
less otherwise specified, we are concerned with the good 
solvent condition in three-dimensional space and adopt 
for simplicity the mean-field value of the critical expo- 
nent v$ = 3/5, where the swelling as against the ideal 
chain 1/3 > 1/2 is due to the excluded-volume interac- 
tions among segments. Increasing the volume fraction 
<p of the segment, one reaches the point </)( ot ') ~ jV -4 / 5 
above which different polymers start to overlap in space. 
In this so-called semidilute regime, the excluded-volume 
interactions are screened at the scale of the correlation 



length £ 



3 / 4 , and the chain conformation can be 
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envisioned as a random walk of the correlation blob of the 
size £ ~ ag 3 / 5 with g segments inside, so that the overall 
size is given by R ~ ^(N/g) 1 / 2 ~ a0~ 1/8 7V 1/2 . In the 
melt, which is a dense limit of the semidilute solution 
0—^1, the excluded-volume interactions are screened at 
the segment scale, so the chain follows the almost ideal 
random walk statistics. 

Turning our attention to the melt of ring polymers, 
the excluded-volume interactions would be screened as 
well, but for this very reason, the topological constraints, 
which may be buried under the excluded-volume effect in 
other situations, stand out. Therefore, unlike the linear 
polymer case, we encounter increased complexity instead 
of simplification due to the screening of the excluded- 
volume effect. Indeed, it has been acknowledged in re- 
cent years that understanding conformations and dynam- 
ics of the melts of noncatenated ring polymers is one of 
the remaining challenges in the polymer science [l0rll2l |. 
Towards this ambitious goal, the construction of a phe- 
nomenological theory would offer promising prospects, 
in which the topological constraints could be accounted 
through the geometrical picture, just as in the tube model 
and the reptation model in the theory of linear polymer 
melts. The aim of the present paper is to provide a de- 
tailed account of the recently proposed such an attempt 
for the statistics of ring polymer melts . 

The rest of the paper is organized as follows. In ScclLTl 
we start with some preliminary remarks on the topologi- 
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cal length scale, and then introduce a mean-held model. 
Sec. IIIII summarizes the results including (i) conforma- 
tional properties of individual rings, (ii) global geomet- 
rical structure of melts and (hi) the effect of the chain 
stiffness along with the comparison with available numer- 
ical/experimental data. We shall also look into previous 
conjectures in light of the present theory. Finally, we 
close in Sec. IIVI with summary and perspectives. 



II. MEAN-FIELD THEORY 
A. Topological length scale 

Entanglements formed in linear polymer melts and 
topological constraints in ring polymer melts differ in 
their nature, but they share some similarities which may 
trace back to the fact that both arise from the non- 
crossability constraint in flexible polymer systems. 

It is well known that in melts of short linear polymers 
TV < TV e , there arc no signs of entanglement, and their 
dynamics is well described by the Rouse model 0-01 • The 
entanglement length N e depends on the chain stiffness, 
concentration, etc., but its typical value in melts of flex- 
ible polymer is on the order of TV e ~ 100 There is a 



quantity similar to TV e for the ring polymer melts, which 
we call the intrinsic topological length and denote by TVi 
in the following. It has been observed in numerical sim- 
ulations that in melts of short ring polymers TV < TVi, 
there seems to be no sign of the topological constraints, 
and each ring obeys the ideal statistics with the size 
R ~ aTV 1 / 2 [l!miii{l2l. There are some good evidence 
that the numerical value of N\ is comparable to that of 
N e in the corresponding linear chain melts [n], [12] • We 
provide a geometrical interpretation of TVi in Sec. IIIIBI 
In linear polymer melts with TV 3> TV e , chains are 
mutually entangled, where the tube size is given by 

1 12 

£ e ~ aN e ' . The tube size and, therefore, TV e as well 
are independent of TV for TV 3> TV e , which are impor- 
tant quantities in the rheological characterization of the 
system. However, the analogy may break down for ring 
polymer melts with TV 3> TVi due to the non-trivial ring 
conformation in large scale in contrast to the ideal con- 
formation in linear polymers. Since there is no reason a 
priori to claim the constant topological length, we let it 
be TV dependent and denote it by g§ (TV) with the "initial 
condition" <7^(TVi) = TVi. The functional form of <?ti(TV) 
and the conformation of individual rings at large scale 
r > alg^N)} 1 / 2 for melts of high molecular weight rings 
TV > TVi should be determined from the topological con- 
straints in the system. This can be done self-consistently 
by the following mean-field theory. 



B. Excluded- volume analogy and Order parameter 

We are concerned with the topological effects in melts 
of noncatenated, unknotted rings. All other effects are 



irrelevant, which leads to the ideal behaviors in melts 
of linear polymers/phantom ring polymers. Two contri- 
butions could be identified. One is the inter-molecular 
(avoiding to link with other rings), and the other is the 
intra- molecular (avoiding the self- knotting) . To write 
down the free energy for these effects, we recall the obser- 
vation that two unlinked ring polymers with zero thick- 
ness, i.e., no excluded- volume, repel each other at close 
distances d, [l9j]. The second virial coefficient of inter- 
action between such rings is found to be of the order 
of the size of the ring cubed. That is to say, in di- 
lute solutions, the topological effect restricting the phase 
space volume can be effectively treated as the excluded- 
volumc problem. There is also a conjecture based on 
the analysis of Gauss linking number of random flexi- 
ble rings that the topological constraints in long flexible 
rings and the excluded-volume interactions have similar 
effects [20j . We shall seek for a possible way to imple- 
ment this excluded-volume analogy in the statistics of 
melts and concentrated solutions. 

In what follows, we include the concentration effect 
by regarding a semidilute solution as a melt composed 
of the correlation blob of size £ ~ a<T> -3 / 4 . The melt 
is understood to be the limit <f> — > 1. The equilibrium 
spatial size R (such as a radius of gyration) of each ring 
is a function of TV and the monomer volume fraction </>. 

Let us introduce several quantities relevant for the 
characterization of ring polymer melts. The total num- 
ber of monomers in the region of volume ~ R 3 of a single 
ring is ~ i? 3 </)/a 3 , thus, the number of rings there is 



TV fi 



a 3 TV 



(1) 



For a given ring, Nr can roughly be thought of as the 
number of neighboring rings, so we call it the coordination 
number, which turns out to be an important measure for 
the geometrical characterization of the melts. Equally 
important is the self-density 



a 3 TV 



(2) 



which measures the degree of compactness of individ- 
ual rings, or the degree of mutual ring intcrpenetration. 
These two quantities are related as Nr ~ (j)/4> s - 

The effective excluded-volume vr of the ring with the 
spatial size R scales with volume and we write it as 



v R ~ R 3 Y 



(3) 



where a dimensionlcss factor Y (independent of TV) is 
introduced to account for the fact that the ring is not 
a rigid particle. Instead, it is a very soft object, and 
significant interpenetrations should occur in melts. Of 
course, this factor Y cannot be arbitrary, and shall be 
fixed later by the physical requirement associated with 
the intrinsic topological length TVi mentioned in Sec. Ill Al 
The "volume fraction" of topological origin then follows 
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as 



R 3 



N R Y 



(4) 



which plays a role as an order parameter in the theory. 
We thus write down the free energy as a function of <fi R 
as 



F{.4>r) = Finter{<t>B.) + Fintrai^B.) 



Noncatenation constraint 



(5) 



Evaluating the noncatenation constraint contribution 
to the free energy is the most difficult part in the 
problem. We follow the excluded-volume analogy and 
make a simple proposal [l3||. Given an effective vol- 
ume fraction <p R of topological origin, we assume that 
its impact can be treated within the framework of liq- 
uid theory, where the excluded-volume effect in dense 
solutions has been intensively studied. We adopt a 
classical van der Waals theory of fluids whose free en- 
ergy density for one component fluid with the vol- 
ume fraction <j> and the excluded-volume v is given by 
f{ft>)/{k B T) = Wv] In [<f>/(l - </>)] - e0 2 with k B T be- 
ing thermal energy. In our case, this may be trans- 
formed as fi„ter{4>R)/{kBT) ~ -[4> R /v R ]]n[{l- 4> R )} 
where we set e = (athermal). The ideal gas term is 
irrelevant here, since it is already included in the non- 
topological part of the free energy. The free energy per 
ring F inter = f inter x R 3 /N R is thus 



k B T 



= - In (1 - <j> R ) 



(6) 



D. Unknotting constraint 



An inspection of eq. ([6]) indicates that increasing the 
ring size R is unfavorable due to the noncatenation con- 
straint. Thus, this topology effect leads to the squeezing 
of the ring toward a globular state, which should be ne- 
gotiated with the unknotting constraint. To derive the 
expression for Fi n t ra {<f>R), it is instructive first to recall 
the free energy cost AFe to confine an isolated linear 
polymer into a small closed cavity under O condition, 
where v$ = 1/2. A standard scaling consideration sug- 
gests the following form; 



D J 



(7) 



where i?e = aN 1 / 2 is the unperturbed coil size and 
D < Rq is a linear dimension of the cavity. The ex- 
ponent Sq can be determined by examining the system 
size dependence of A.Fq under constant segment density 
(or constant osmotic pressure). That is to say, we re- 
quire A_Fe — > kAF& upon enlarging the system size as 



N — > kN and V ~ D 3 — > kV. This thermodynamic 
consideration fixes Sq = 6 H, [2l[ . 

We expect that the probability of self-knotting depends 
on the self-density <p s , and it linearly increases with the 
system size N under constant <f> s . Therefore, the above 
reasoning for AFq can be directly applied to Fi n t ra , too, 
so we obtain 



F l : 



k B T 



R J \RJ 



3 ^>- 3 / 4 



N 4> 



(8) 



where Rq ~ a^ -1 / 8 ^ 1 / 2 is the unperturbed polymer size 
in semidilute regime as already stated in Sec. HI and R < 
Rq is a linear dimension of individual rings squeezed by 
surrounding rings. Using the definition of 4> R (eq. (|^])) 
and N R (eq. ([T])), the above equation can be rewritten as 



F intra {<t>R);N NY 2 cf> 5 / 4 



k R T 



(9) 



which guarantees the linear dependence on N under the 
fixed order parameter <j) R . 

To get a physically appealing insight, let us again come 
back to the problem of the confined O-chain. Its spatial 
structure can be viewed as a dense piling of blobs of size 

1/2 

£e — a 5e with ge segments inside. This implies the 
relation a 3 ge/^e — (f> — a 3 N/R 3 4^ £q ~ a^ 1 and 
gs ~ cj)~ 2 . The free energy cost (eq. (O) can be obtained 
by assigning ~ k R T per blob, i.e., AFs/ksT ~ N/gQ ~ 

It is tempting to follow the same line in our problem, 
too, but now the self-density cf> s should be used to write 
down the space-filling condition. Let us do so. A view 
of the compact ring of size R (< Rq) as a dense piling 
of topological mesh (blob) implies ^ ~ ^(g^/g) 1 ^ 2 and 
a3 g^/^ = 4>s- One thus finds 



<h 



i^ 3/4 N R ~ CN R 



3/4 s - 



r b/i N 2 R 



gN 2 R . 



(10) 

(ii) 



Assigning 

F in tra/(k B T) ~ N/g^ ~ R 3 /^, which coincides with 



k B T per topological strand, we obtain 



intra 

the scaling derivation of F intr a (eq. ©). Moreover, the 
deduced relations ^ ~ £N R and g§ ~ gN R point to an in- 
teresting relation between the topological constraint and 
the geometrical property of the melt. Indeed, it indicates 
that the number of units (correlation blobs) required to 
form a topological constraint is of order ~ N R . 



E. Short ring case and Y-factor 

To complete the formulation, the factor Y introduced 
in eq. ([3]) remains to be fixed. This can be done by con- 
sidering the behavior of the short ring melts with N < N% 
(or the short scale behaviors of longer ring melts). 

To begin with, let us recall the logical structure of the 
mean-field theory. The free energy given in eq. ([5]) is a 
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so-called constraint free energy obtained by partial sum- 
mation over phase space with <pR (order parameter) fixed. 
The true thermodynamic free energy T is given by 



exp [— (3 J 7 } 



b R exp[-/3F(fo)] 



where (3 = (kgT) 1 . For long ring solutions (N ^> Ni), 
i.e., large systems, the above integral is dominated by the 
narrow region around the minimum <J)r 
to 



n q , leading 



and 



4 e9) 



represents the 



from which N 



equilibrium value of 

£*> = ^/Y and ~ 

aiV 1 / 3 (^ e9) ) 1 / 3 ^- 1 / 3 follow; see eq. © and ©. 
Note that superscript (eg) marking the equilibrium 
value will be attached only when it is necessary, and 
will be omitted unless confusion is expected. However, 
for the melts with short rings, i.e., small systems, 
fluctuations around the minimum arc substantial, and 
the above minimization procedure is insufficient to 
capture the equilibrium properties of the system. 




FIG. 1: (Color online) Free energy profiles as a function of 
4>r for the melt with short rings N/g* = (0.4) 2 (left) and long 
rings N/g* = 10 (right). 

For example, we plot in Fig. Q] (left) the free energy 
profiles for the melt of short rings (N ~ iVi), where the 
free energy minimum is attained at low (f>R < F an d the 
corresponding free energy is comparable to the thermal 
energy. Physically, this means that for melts with shorter 
rings (JV < N\), the topological constrains are buried in 
the thermal fluctuations, therefore, ideal chain behavior 
is expected there. One may find a situation Fi nter (4>R) ~ 
Fintra{4>R] N\) ~ fcfiT at (f>R. From this, we find the 
expression for the factor Y; 



Y ~N- 



-1/2,-5/8 



Nl 

9 



-1/2 



f!2) 



where <f>R is an onset value of the order parameter above 
which the topological constraints set in. The value of 4>r, 
which must be smaller than unity, is a universal number 



in our theory in the sense that its precise value is in- 
sensitive to the structural details of individual polymer 
systems, rather determined by a geometrical requirement 
(see Sec. IIII B[) . Figure [1] (left) indicates a reasonable 
range for it around (/>r ~ 0.5. Later, in Sec. IIIIB1 we 
shall see that 4>r — 0-4 can indeed account for recent 
numerical simulations very well. In the last expression, 
£i ~ ii.Nx/g) 1 / 2 ~ aNl ,2 (j)- 1 ' s ' is the spatial length scale 
corresponding to N\. 

As expected, the Y-factor depends on the intrinsic 
topological length Ni and acts as a softening factor com- 
pared to the rigid particle, which allows the partial in- 
terpenetration of different rings. By lumping Ni and 4>r 
together, one can rewrite cq. (|T2|) as 



where 



Y 



<■* _ £i 

4>R 



-1/2 



1/2 



3 =12" 



(13) 



(14) 



are the size of the topological mesh and the number 
of segments in it, respectively, in the compact statistic 
regime as detailed below in Sec. IIIIl 

The intra-molecular free energy eq. ^ can now be 
rewritten as 



intra 



k B T 



N 



(15) 



III. RESULTS AND DISCUSSION 

Now the stage is set, so we proceed to seek for the 
spatial structure of rings in large scales N > Ni. In the 
following, whenever concrete numerical calculations are 
necessary, we set all numerical prefactors of order unity 
(indicated by the ~ symbols in equations) to be unity. 
As shown in Fig. [T] (right), here, the energy scale is much 
larger than the thermal energy, thus the equilibrium size 
can readily be determined by minimizing the free energy 
with respect to </>#. 



A. Conformation of individual rings 

In Fig. [21 we plot the equilibrium size R of rings as a 
function of the polymerization index N. The following 
trends can be read off; (i) sufficiently long rings (N > N2) 
obey compact statistics with v = 1/3, (ii) rings with in- 
termediate length {N\ < N < N2) are well characterized 
by the effective exponent v e ff — 2/5. Therefore, one 
can say that the statistics of ring polymer melts have 
two well-defined regimes: ideal and compact statistics 
for short N < Ni and long N > N 2 rings, respectively, 
which are separated by a wide crossover region spanning 
over more than one order of magnitude. 
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N/g* 



FIG. 2: (Color online) Master curve of the ring size R/£* as a 
function of the chain length N/g* (double logarithmic scale) 
obtained by numerically minimizing the free energy eq. ([5]). 



To get an insight into this crossover, we note that from 
eqs. (HJ, flD, @, (fT2 ]) .(fI5 ]l . a mean-field solution for TV > 
Ni can generally be written in the form 



where £* = a 



R = C 



= «rA- 1 /8 £ ,*i/2 



A^ 



1/3 



(16) 



60a 1 (eq. d), and de 



notes the equilibrium value of the order parameter which 
depends solely on N/g*, i.e., 4>^ q) = ^ (N/g*). The 
above equation can be transcribed into a physically ap- 
pearing form as 



R = ^(N) 



N 



1/3 



(17) 



on N is absorbed 



where the subtle dependence of <^ R 

in the topological length as £tj(A) = £*(j) R 3q> and 9h{N) 

g*(4>ji ) 2 - This indicates that individual rings arc 
squeezed by the noncatenation constraints, and hence, 
follow the compact statistics at fixed N. In addition, 
unlike the usual globule found in liner polymers, the un- 
knotting constraint of rings leads to a unique conforma- 
tion known as a crumpled globule (24|, i.e., upon renor- 
malization a large scale spatial organization of individ- 
ual rings is formally characterized by a space-filling curve 
with fractal dimension Df = v^ 1 = 3 (see Sec. lIIIE|) . As 
this expression shows, the peculiarity in the ring con- 
formation in melts lies in the fact that a renormaliza- 
tion factor, which may be called topological mesh [24| . 

£tj(iV) = ag^ 2 (N)cj)^ 1 ^ is length dependent as a result 
of the competition between inter-molecular and intra- 
molecular topological constraints. 

Looking at the behavior of the order parameter, it 
takes the value 0^ 



constraints (N = N\), and increases with N towards the 

saturation — > 1 at N — > A 2 , where the saturation 
point N = N2 can be identified as the onset of the com- 
pact statistics. Correspondingly, g§ and ^ change in such 
a way that = N\ and ^ = £1 at N = Ni and they in- 
crease with N towards the saturation g§ = g* and ^ = £* 
at N = A 2 . As N 2 > g* = Ai</>^. 2 (a formula to estimate 
A 2 will be given in eqs. (j2"5j) . (fMf and (j3"rjj) along with 
its physical interpretation), the evolution of g^, ^ during 
this interval is quite slow. Indeed, this slow evolution is 
controlled by the logarithmic divergence towards a topo- 
logical dense-packed limit — > 1), but well fitted by 
the effective (Active) power laws ^ ~ A 1 / 5 and ^ - A 2 / 5 
(Fig. [3]). This renders a small correction to the exponent 
v, which leads to the effective exponent in the ring size 
R ~ <^(A/c^) 1/3 ~ A 2 / 5 . Such a feature has been ob- 
served in several numerical simulations [HI, [IH, [25l - [27| 
and real experiments (38j with moderate chain length, 
which indeed coincides with the Catcs-Dcutsch conjec- 
ture Q. However, we emphasize that v e ff ~ 2/5 found 
here is not a well-defined exponent in our theory, but a 
fictive exponent seen in a broad (but restricted) range of 
the chain length Ai < A < A 2 . 




bn at the onset of the topological 
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FIG. 3: (Color online) Topological mesh size as a func- 

tion of the chain length N/g* . (Left) is a semi-log plot show- 
ing a logarithmic increase towards to ^ — > £*, while (Right) 
is a double logarithmic plot showing an effective power low 
£ q ~ N 1 ^ 5 in the intermediate range of the chain length 
(N\ < N < N 2 ). 

For large enough molecular weight (A > A 2 ), the 
compact statistics R = ^(N/g*) 1 ' 3 is realized in accor- 
dance with recent numerical and experimental observa- 
tions [ll|, [TH, [I?], [U| ■ The coordination number follows 
from eq. flU as N R ~ g* 1/2 5 / 8 , and thus the self-density 
4> s ^ <t>/N R ~ g*' 1 / 2 ^/* (eq. ©). This agrees with re- 
sults for the concentrated solution of randomly branched 
polymers derived by Daoud and Joanny (28l . |29[ | . Their 
expressions for the chain size, the coordination number 
and the self-density in melt limit 0—5-1 are, respectively, 
R ~ a N^ d A^- d y 2d , N R ~ A( 2 - d )/ 2 and S ~ A^ 2 )/ 2 
for the space dimension 2 < d < 4, where A 2 is called 
the activity of the trifunctional units (branching points) . 
Noting that A -1 ~ g* corresponds to the linear parts 
between trifunctional units, our results are recovered by 
inserting d = 3. It is interesting to observe that Nr 
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(the degree of the overlap) is controlled by A (the degree 
of the branching). Indeed, Nr increases with decreas- 
ing A, and crosses over to the value of linear chain melts 
N R ~ AfW 2 )- 1 a t A ~ AT -1 , where a labeled chain is sur- 
rounded by a large number of other chains, which grows 
as iV 1 / 2 in d = 3, leading to significant entanglements. 

The optimum free energy in this compact statistic 
regime is F/(k R T) ~ N/g§ , i.e., an order of the mag- 
nitude of thermal energy per topological strand. This 
translates into the osmotic pressure n toJ) ~ FN R /R 3 ~ 
k B T<j>/(d?gfi ~ (N R )" 2 k B T(f, 9 ^/a 3 . Adding this to the 
usual excluded-volume contribution IIo ~ k R T/^ 3 ~ 
fcsT(/) 9 / 4 /a 3 0, the total osmotic pressure thus can be 
written as 



n = no + n top ~(i + A^ 2 )n . 



(18) 



It would be instructive to compare this result with the 
change in the osmotic pressure in linear chain solutions 
due to an increment in the local excluded-volume. Let 
v ~ a 3 and c be the excluded-volume and the concen- 
tration, respectively, so that <j> = cv. Then, the os- 
motic pressure is Ho/k R T ~ v -1 ^ 9 / 4 = ii 5 / 4 c 9 / 4 . Let 
us slightly increase the excluded-volume as v — » v + 5v 
(dv/v « 1). The resultant increase in osmotic pressure 
is SU (v) = Ho(v + 6v) — n (u) ~ TLq(v)8v/v. Therefore, 
the effect of topological constraints seems to be reflected 
as an increase of the local excluded-volume interactions, 
i.e., N R 2 O Sv/v in accordance with the des Cloizeaux 
conjecture [20| . By inserting the value of the coordina- 
tion number estimated below (Sec. IIII B[) . we see that 
the additive term is too small to be detected in practical 
measurements 121. 



B. Large-scale spatial structure 

During the broad crossover region (A r i < N < N2), the 
melt of rings adjust its own spatial structure by balanc- 
ing the noncatenation and unknotting topological con- 
straints. The coordination number N R (eq. (^Q)) is a good 
measure to capture the large scale (r > R) geometrical 
feature of the melts. According to eqs. (U]) and (fT2} . it is 
expressed as 



N R (N) = 



Y 



Ni 



1/2 



(19) 



so that it slowly increases from N R (Ni) at N — Ni to- 
wards N R (N 2 ) at N = N 2 . These two values N R (Ni) 
and N R (N 2 ) turn out to be key quantities in the follow- 
ing discussion. 

• At the onset of topological constraints (N = Ni): 
We have c/>r(Ni) = <fi R so that 



Nj 
9 



1/2 



1/2 ,5/8 



(20) 



Hence, Nr(Ni) is found to be directly linked to the 
intrinsic topological length N±. A rough estimate 
Ni ~ 100, 4> ~ 0.5 gives N R {N X ) ~ 6. 

• At the onset of compact statistics (N = N2): We 
have 4>r(N2) = 1 so that 



Nr(N 2 ) = I 



N R (N! 



1/2 



g *l/2^/8 (21) 



A rough estimate Nr(Ni) ~ 6, 4>r ~ 0.4 gives 
Nr(N%) ~ 15, which is in reasonable agreement 
with recent numerical simulations (see the caption 
of Fig. [4}. Its geometrical interpretation was sug- 
gested in ref. [11|, where the correspondence of 
Nr(N 2 ) with the coordination number of jammed 
hard ellipsoid systems [23| was pointed out. Hence, 
Nr(N 2 ) is supposed to be linked to the geomet- 
rical property inherent in the packing in three- 
dimensional space. 

According to cq. (f2"Tj) . the above consideration identi- 
fies the Y- factor ~ 1/N R (N 2 ) to be a geometrical quan- 
tity, too, which in turn suggests the geometrical origin of 
the intrinsic topological length scale N\ (eq. (fT2")l ). thus 
also N R (Nx) as well as g* (cq. dTgl)). 
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FIG. 4: (Color online) Dependence of the coordination num- 
ber Nr on the chain length N. The solid line is the the- 
oretical prediction, while points are results from numerical 
simulations (read off from Fig. 8 in ref. p4| . and Fig. 13 in 
ref. P3]). Circle: Monte Carlo simulation by Vettorel et. al. 
(N e ~ 175, N R (N 2 ) ~ 15) [H]], Square: Molecular dynamic 
simulation by Halverson et. al. (N e ~ 28, N R (N 2 ) ~ 17) [12]. 

In Fig. we compare the prediction of our mean- 
field theory with the results from recent numerical sim- 
ulations. To determine the scaling factor g* , we first 
identify the intrinsic topological length N± to be equiva- 
lent to the entanglement length N e which was measured 
in simulations, and then assigned 4>r = 0.38 for the or- 
der parameter at which the topological constraints set 
in. This yields g* 
g* ~ 190 for 



1200 for the model in ref. 
that in ref. 



11] and 



12J. It is remarkable that 
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with the above procedure for the rescaling, i.e., (i) the 
identification of Ni — A e , (ii) the assignment of a rea- 
sonable and common value of 4>r, two different sets of 
simulation data are almost perfectly collapsed into the 
theoretical master curve. The fact that the common <j)R 
value can fit the results from different simulation mod- 
els (flexible rings in lattice Monte Carlo simulation 11 1 
and semi-flexible rings in off-lattice Molecular dynamics 
simulation [l2|) indicates that 4>r is insensitive to the 
model details. This is an expected result from the fol- 
lowing consideration; From eq. (j4} with <p^' {N\) = 4>r 



and 4> [ n q> (N 2 ) = 1, it follows that 

t _ N R {N{) 
PR N R {N 2 ) 



(22) 



As already described above, Nr{N\) and Nr{N 2 ) are ge- 
ometrical quantities inherent to three-dimensional space, 
so their ratio should be universal, too. 



C. Physical picture of the crossover 

After clarifying both the intra- and inter-molecular 
spatial structures, we now ask what is the underlying 
physical picture of the broad crossover regime (N± < 
A < N 2 ). Let us look at one given ring (shaded circle 
in Fig. [5]) in the melt and call it a reference ring. Define 
the number 



m(A) 



N 



(23) 



of the topological meshes per ring. It takes m(Ni) — 
Ni/g^(Ni) = 1 at the onset of the topological effect 
N = Ni, and slowly increases to m(N 2 ) = N 2 /g^{N 2 ) = 
A2/.9* at A = N 2 . At even larger A > A 2 , g^(N) satu- 
rates, so that m(N) linearly increase with A in the com- 
pact statistic regime. These meshes are forced to form 
due to the intra-molecular topological repulsion against 
the squeezing by surrounding rings. Therefore, any topo- 
logical mesh in the reference ring is under the interference 
of the surrounding rings, as schematically illustrated in 
the left column in Fig. [5] 

Let us further define the occupation rate 



M(N) 



N 



m(A) 

= N R (N) ~ g[N R (N)]f 



(24) 



where the last relation follows from eqs. (fTTj) and (|23| . 
It represents the fraction of the Nr surrounding rings 
which take part in the interference with the reference 
ring, thus forming an active pair with the reference ring 
in the sense of the topological constraint. At N = Ni, 
we have M{N X ) ~ l/N R {N x ) < 1. This means that at 
the onset of the topological effect, there is some freedom 
to form a nascent mesh, i.e., it can be formed coopera- 
tively, where all N R (Ni) rings moderately contribute to 
the net interference energy ~ fcgT, or alternatively it 



can be formed by an intense interaction (on the order of 
~ ksT) with any one of the N R (Ni) rings. In this situ- 
ation, there is some combinatorial freedom to create the 
mesh. 

Now suppose the hypothetical process in which we in- 
crease the chain length A with <f> R fixed. It is clear from 
the free energy eqs. ©, (fT5)) with their behaviors shown 
in Fig[T]that the created state is unstable with high intra- 
molecular energy cost F intra with lots of meshes. Switch- 
ing off the constant 4> R constraint, therefore, the system 
relaxes to the equilibrium state with larger <f> R by enlarg- 
ing the ring size R at the expense of the intcr-molccular 
energy cost Fi nter with the larger coordination number 
Nr. The size of meshes ^ also grows with N so that 
the self-density <j> s ~ a 3 g^/£ 3 ~ { a /^\\)4> 1 ^ decreases. 
This dilution of the self-density is caused by the inter- 
penetration of the surrounding rings, i.e., the increase 
of Nr. In essence, the topological constraints become 
tighter with increasing A as a result of the self-consistent 
spatial structural adjustment. 

In fact, by comparing the growth rate of m(N) and 
Nr(N) (adopting Active exponents obtained in Sec. lIII Al 
these quantities scale as m (A) - A 3 / 5 and N R (N) ~ 
A 1 / 5 in the crossover regime), one finds that the occu- 
pation rate M(N) ~ N/(g[N R (N)} 3 ) grows with A, and 
eventually reaches the saturation point, i.e., the onset of 
the compact statistic, which can be defined by the rela- 
tion M(N 2 ) — 1. From this condition, we find 



(25) 



where the last relation follows from eq. (|20[) . Substituting 
the estimates N R {Nx) ~ 6, N R {N 2 ) ~ 15 and Ai = A e 
from the simulation data, we find the estimates N 2 ~ 10 4 
for the model in ref. [Tl| and A 2 - 10 3 for that in ref. [13], 
in rather good agreement with observations. 

Physically, this saturation condition means that there 
are Nr meshes in the reference ring, all of which inter- 
fere with one of the surrounding Ar rings. There is, on 
average, one-to-one correspondence between mutually in- 
terfering meshes. Hence, the crossover regime is charac- 
terized by the continuous loss of the combinatorial free- 
dom to form the network of meshes, and M(N) ~ 1 at 
A = N 2 signifies the densely packed (no voids) state in 
terms of the topological volume fraction 4>r — > 1. 

Figure [5] summarizes this subsection. At A > Ai, 
the system of ring polymer melts can be viewed as com- 
posing of topologically active clusters. Rings belonging 
to the same cluster cooperatively create meshes within 
these rings. If we look at one ring (reference ring marked 
in the middle column) in a particular cluster, some of its 
neighboring rings belong to the same cluster, but oth- 
ers do not. The increase of the occupation rate M(A), 
which is the fraction of the former, with A corresponds 
to the growth of the cluster size. At M(N = N 2 ) ~ 1, 
the combinatorial freedom for the clustering patterns is 
lost, and the whole system is spanned by a single cluster. 
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R(N)=^(N) 



[ gl <N<g 3 ] 
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[N=g 3 ] 
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FIG. 5: (Color online) Schematics illustrating the spatial 
structures in ring polymer melts, [top] almost ideal ring 
regime (N < Ni), [2nd line] onset of the topological con- 
straints (N = ATi), [3rd line] crossover regime (Ni < N < 
N2), [bottom] onset of the compact statistics (N = A/2)- 
The right column shows a reference ring (shaded) and its 
neighbors, which are represented by soft spheres. During the 
crossover, the coordination number Nr, thus, the degree of 
the overlap measured by cj>n increases with N as shown in 
the right column. The middle column represents the global 
structure focusing on how the size and the correlation of topo- 
logically active clusters grow with N. The left column shows 
the internal structure (meshes represented by dotted circles) 
of individual rings, and the finer inter-mesh structure with 
(a part of) the reference ring (thick line) and a part of other 
neighboring ring (thin line) penetrating the reference ring re- 
gion. 



ring overlap . We shall examine such a stiffness effect 
here. 

With the increase in the local chain stiffness, each sta- 
tistically independent segment becomes anisotropic with 
the thickness a and the length I > a. We define the poly- 
merization index as N — L/a, where L is the contour 
length of the ring. The number of statistical (Kuhn) 
segments per ring is thus Nk = L/l = N/p, where 
p = l/a is the stiffness parameter (aspect ratio of the 
segment). According to the standard theory of polymer 
solutions, there are two qualitatively different regimes 
which are separated at the threshold stiffness parame- 
ter p* ~ (j)- 1 ! 3 HUGH. At p < p*, the correlation 
along the chain can be treated through the blob as in the 
flexible chain solutions, which has now p dependence as 



3/4 p- 1/4 , 
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-5/4-3/4 



(p<p*)- (26) 



The unperturbed (i.e., no topological constraint) polymer 
size which appears in eq. (HJ) is thus i?o — ^(N/g) 1 ^ 2 ~ 
a 0-i/8jvi/2pi/8^ On the other hand, &tp > p*, the corre- 
lation along the chain is negligible and the mean-field the- 
ory is reliable there. In this mean-filed regime, the chains 
without topological constraint exhibit nearly ideal statis- 
tics on all length scales larger than the segment length I, 
where it is convenient to define the quantities 



P 



aG 1 ' 2 p 1 ' 2 1 G~0" 



V 3 



(p>p*W) 



so that R ~ -(iV/G) 1 / 2 ~ aiV 1 / 2 ^ 1 / 2 . In subsequent 
formulas below, these quantities play analogous roles as 
g and £, respectively, in this mean- field regime p > p*. 

Following the construction of free energy in Sec. Ill CI 
and HID I it is easy to see that eq. ([6]) for Fi nter remains 
intact, but the appearance of eq. (JSJ) for Fi ntra is modified 
due to the p dependence in Rq] 



J 1 ^ ip>P) 



By following the argument in Sec. Ill El we can fix the Y- 
factor from the condition F mter (4> R ) ~ F intra ((t>B,; Ni) ~ 
k B T; 

y ~ { ^r v V- s/8 P- 3 / 8 = (f)" 1/a ^ (p<p*\ 29) 

4 ~^ R N^ 2 <j>- l p-^={W 1/2 ^ (P>P*) ' 



which can be rewritten in the form of eq. (|13[) as 



D. Effect of stiffness 

In linear chain melts, the chain stiffness is known to 
be an important factor controlling the degree of entan- 
glements, i.e., a slight increase in the stiffness promotes 
it by reducing N e substantially [l4|. Similarly, in the 
ring polymer melts, it is known that tuning the persis- 
tence length provides an efficient route to increase the 



Y 



with 



6_ 

'PR 



,s -1/2 
B-) (P<P*) 
! - 1/2 (30) 
^) (P>P*) 



,\l/2 
9 T ) (P<P*) 
1/2 

(p > p*) 



(31) 
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Substituting this result into eq. (|28j) . it turns out that 
Fintra can be generally written in the form of eq. (fT5|) . 
Therefore, an important p dependence can be absorbed 
in the intrinsic topological length iVi (or equivalently <?*). 
Hence, all the results (including Figs [2j [3j [4j expressed in 
the rescaled forms already encompass the stiffness effect. 

For moderately long ring systems (N > Ni), the topo- 
logical meshes are created. From the space filling condi- 
tion of these a 3 ^/^ 3 ~ 4> s (Sec. Ill I)]) with the intra-mesh 

statistics ^ - £(V.9) 1/2 - a^ 1/8 5 t 1/2 p 1/8 for V < V* 
and £t| ~ E(g ii /G) 1 / 2 ~ ag\j 2 p x l 2 for p > p* , we ob- 
tain the mesh size g^ and ^ in terms of the coordination 
number as follows; 



en r 



(p < p*) 
(p > p*) 



/ gN R ( P <p*) 
\GN 2 (p>p*) 



(32) 



(33) 



which are generalizations of eqs. (JTUJ) and (fTT|) . 

From the condition M(N 2 ) — 1, we can generalize 
eq. (HH]) as 



(34) 



Using the mesh size relation of eq. (|3"3")l at N = N\, i.e., 
g^{N\) = N\ (cf. eq. (f2"U)) ). one can explicitly write down 
the p dependence of N\ as 



TV! ~ [TVflXiV!)] 2 X 1 



(p<p*) 

G ( P >p*) 



(35) 



Then, eq. (|34|) for N 2 can be collectively written as 



N 2 ~ ATi x 



[Afo(JV 2 )] J 



N^ R 2 N R {N 2 



(36) 



which is nothing but the last expression in eq. (|25|) . i.e., 
again the p dependence is absorbed in N\ . 

In Fig. |51 we plot N\ as a function of p according to 
eq. (|35p . together with the data of 7V e obtained in linear 
chain melts [l4[. It is apparent that these two quanti- 
ties share the same qualitative trend. Such a coincidence 
provides a further support for the conjecture iVi ~ N e , 
pointing out an interesting link between the topological 
constraints in ring polymer systems and the entangle- 
ment concept in linear chain systems. 

The expression in eq. (|36[) indicates that the ra- 
tio N2/N1 is a model-independent universal number ~ 
N R (N 2 ) 3 /N R (N 1 ) 2 - 50-100 independent of p. In melts 
and concentrated solutions, the threshold stiffness p* is 
quite small, hence, both onset lengths of the topologi- 
cal effect (TVi) and the compact statistics (iV 2 ) decrease 
sharply with the increase in p. As a consequence, the re- 
gion of the crossover regime N G (Ni,N 2 ) becomes nar- 
rower and the compact statistic regime is reached with 
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FIG. 6: (Color online) Dependence of intrinsic topological 
length scale iVi on the stiffness parameter p. Curves are 
drawn according to eq. I|35p . where numerical coefficients are 
appropriately adjusted so that the two curves coincide at 
p — p* ~ 2. Also plotted is the entanglement length N e 
in linear chain melts determined by primiti ve p ath analysis 
(points and error bars: from table [1] in ref. (141]). 



shorter rings for a semiflexible ring in agreement with 
numerical observations [l2|, EH [H| • Further increase in 
p may lead to the occurrence of the local orientational 
ordering p^j . 



E. Comparison with other models 

Here we shall look into other models of ring polymer 
melts in the light of the present theory. 

Cates-Deutsch conjecture: In their seminal paper, 
Cates and Deutsch (C-D) argued that unconcatenated 
rings in the melt may have statistics intermediate be- 
tween those of compact (y = 1/3) and ideal [y = 1/2) 
chains Q ■ Specifically, they proposed a conjecture on the 
scaling exponent v = 2/5 by assuming and minimizing 
the following free energy: 



F(R) 
k R T 



a 3 N 



Na 2 
R 2 



(37) 



Their reasoning for the above free energy is as follows. 
The first term is the coordination number N R (with the 
melt limit <f> = 1), and they assumed that roughly one de- 
gree of freedom is lost for each of the N R neighbors which 
the ring is prevented from threading, i.e., the noncatena- 
tion constraint. The second term, on the other hand, rep- 
resents the entropic penalty of squeezing an ideal chain 
into a small space of size R < aN 1 ^ 2 . We now argue 
that (i) the first term (inter-molecular effect) in eq. (|3"T)) 
corresponds to the "dilute" limit in our theory, and (ii) 
the second term (intra-molecular effect) is qualitatively 
incorrect. 

(i) First, eq. {[SJ) for Fi nter can be expanded in a virial 
series Fi nter /k R T ~ <fi R + 4> 2 R /2 + ■ ■ ■ . Retaining only 
a lowest order term (second virial approximation), one 
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finds Finter/ksT ~ R 3 (f>Y/ (a 3 N) , which corresponds to 
the free energy adopted in the C-D theory aside from 
a factor Y. Indeed, Catcs and Deutsch were careful 
enough to make comment that the estimated free en- 
ergy Fi n t e r / '(fcflT) ~ Nr due to the nonconcatenation 
constraint may be generally written in the form ~ Nr 
with an unknown exponent a. Our proposal utilizing the 
van der Waals theory provides a possible (yet nontrivial) 
way to gather up all the terms, leading to Fi nter eq. (|6|). 
In addition, it turned out that the higher order terms 
neglected in C-D theory should be important. 

(ii) To see the second point, let us remember the dis- 
cussion in Sec. Ill Dl Since any ring in melts is squeezed 
by surrounding rings from all directions, the situation 
is reminiscent to the chain confined in a closed cavity. 
As we have already seen, the scaling form in eq. (|37|l is 
only valid for an ideal chain. To make the point clear, 
it is convenient to think about a reference linear chain 
in the linear chain melt. If, by some mean, we pres- 
surise this reference chain into a small size R < aN 1 / 2 , 
the associated free energy is indeed given by the sec- 
ond term in eq. (|37|) . However, this does not apply to 
the ring polymer counterpart because of the topologi- 
cal reason. Indeed, the squeezed ring polymer takes a 
special entanglement-free conformation called crumpled 
globule [H| (see below), whose statistical properties are 
rather different from those of a squeezed ideal chain. This 
is manifested, for example, in the profile of the segment 
concentration along the cavity axis for the squeezed ring. 
While the sinusoidal type function describes that of the 
ideal linear chain (or ideal phantom ring) 0], our blob 
picture in Sec. Ill Dl suggests that the ring polymer de- 
velops a bulk plateau profile except for the boundary re- 
gion of the size ^ (eq. (fTUj) ). Such a feature has indeed 
been observed in numerical simulation [22| . Therefore, 
the main free energy penalty associated with squeezing 
the ring polymer arises from the unknotting constraint, 
and this is evaluated as eq. (JHJ, ([9]) or (fTB"j) as we have ar- 
gued in Sec. Ill Dl In summary, both the inter-molecular 
and intra-molecular free energies are underestimated in 
C-D theory, and the exponent v = 2/5 results from the 
cancellation of these two errors. 

Lattice animal model: The idea of reptation in linear 
chain dynamics was first formulated for a single chain in 
the fixed obstacles (cross-linked gel) [f| , and then applied 
to the rheology of entangle polymer solutions Q. There 
have been attempts to follow the same procedure in the 
ring polymer case H [1(3, Ull HH . A long ring polymer 
placed in the array of fixed obstacles (without making 
any concatenation with them) takes a particular confor- 
mation known as the lattice animal whose property is 
well understood [H, [32|, . Indeed, as already stated in 
Sec. IIII A| our results for the compact statistics regime 
(N > N2) coincides with those obtained for the concen- 
trated solutions of randomly branched polymers, which 
are continuous analogues of lattice animals. 

In melts of ring polymers, there are no fixed obsta- 
cles, but each ring feels topological constraints due to 



the presence of surrounding rings, and one may say that 
this could be mapped onto the single ring problem in the 
"mean-field" array of fixed obstacles with their spacing 
^ created by other rings. There is one important missing 
ingredient, however, which should be properly taken into 
account. Although tempting, the topological effect in the 
above picture is undertaken only in the inter-molecular 
interactions. In reality, however, the "mean-field" can 
self-adjust by weakening the inter-molecular penalty at 
the expense of the intra-molecular one (see discussion in 
Sec. IIII Cj) . The present theory provides a physical picture 
of how this self-adjustment is achieved, indicating that 
the topological meshes in individual rings build up along 
with increasing contact number Nr. Hence, the "mean- 
field" obstacles should be assumed to originate from both 
inter- and intra-molecular topological constraints, and 
the resultant conformation of individual rings is remi- 
niscent to that of crumpled globules [24| (see below). 

Crumpled globule model: The concept was originally 
proposed by Grosberg, et. al. as a long-lived kinetic 
intermediate on the pathway of the linear chain col- 
lapse [23 |. and then, hypothesized as a large scale DNA 
organization in interphase chromosomes |36[. More re- 
cently, its relevancy to the conformation of individual 
rings in their melts has also been discussed flU [l2j . 

To see the point, let us revisit our discussion in 
Sec. Ill Dl where we derived the free energy term Fi ntra 
due to intra-molecular topological constraints in analogy 
with the confined polymer problem. When confined in a 
small closed cavity with the size D smaller than the bulk 
coil size Rq ~ aN V3 , the polymer takes a compact globule 
conformation. For a linear polymer, it can be viewed as a 
piece of melts, or more precisely scmidilute solutions with 
the mesh size £ ~ a(a/ £>)3"3/(i-3* 3 ) N u 3 /(i-3„ 3 ) [7j [2ijj 

It implies that there are many intra-molecular entangle- 
ments (self-knotting), and the Flory theorem tells that 
the conformation of the sub-chain is characterized by the 
trajectory with the fractal dimension Df = 2 in the scale 
£ < r < D because of the screening of excluded-volumc 
interactions. However, such a self-knotting is forbidden 
for ring polymers, and upon confinement in a cavity, it 
follows a self-similar space- filling trajectory with Df = 3, 
that is the crumpled globule. Here, all chain parts on 
scales larger than the topological mesh ^ are mutually 
segregated in space due to the topological constraint. 

Our intra-molecular free energy Fi ntr a indeed follows 
from the crumpled globule concept, which combined with 
the inter-molecular squeezing effect Fi nter naturally leads 
to the conformational properties discussed in Sec. IIII Al 
But again, there is one important difference compared to 
the above example of rings under external confinement, 
where the topological mesh size is assumed to be constant 
(usually assumed to be ^ ~ £,(N e / g) 1 / 2 ). The same ap- 
plies to the problem of the concentration solutions of ran- 
domly branched polymers [28[ (see Sec. IIII A[) . in which 
the corresponding scale ^ is set by the degree of the 
branching A, that is a quenched variable. In contrast, 
as explicitly written in eq. (|17|) . the topological mesh 
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size ^ in the melts problem is scale dependent, which 
should be self-consistently determined by the balance be- 
tween Fi n ter and Fi n tra- This leads to a broad crossover 
regime with the Active exponent v e ff — 2/5 during which 
the topological constraints of the system becomes tighter 
and tighter, as we have already seen in Sec. MI Al IIII Bl 
andETCj 



IV. SUMMARY 

In this paper, we have sought for a geometrical repre- 
sentation of topological effects in melts and concentrated 
solutions of ring polymers. Our central idea lies in the 
introduction of the topological volume fraction, through 
which the noncatenation topological constraint is imple- 
mented via the effective excluded- volume effect. From 
this, we have constructed a mean-held theory, which pro- 
vides us with a self-consistent topological constraint de- 
termined by the balance between inter-molecular non- 
catenation and intra-molecular unknotting requirements. 

From the beginning, we have assumed the presence of 
an intrinsic topological length scale N\, which signals 
the onset of the topological effect, thus playing a role 
analogous to the entanglement length N e in the rheol- 
ogy of linear chain melts. Once N\ is fixed, our theory 
identifies two additional length scales g* = Nicj)^ 2 and 
N 2 = N 1 4>^ 2 N R {N 2 ) = g*N R (N 2 ). The latter signals 
the onset of the compact statistic regime, while the for- 
mer represents the number of segments in the topological 
mesh in that regime. 

The geometrical interpretation of Ni is provided via 
the coordination number N R (N). We have seen that its 
values N R (Ni) at TV = JVi and N R (N 2 ) at N = N 2 are in- 
deed geometrical quantities associated with the random 
packing properties in three-dimensional space. There- 
fore, the ratios g*/N x = 4> R 2 = [A^A^/A^Ni)] 2 and 
N 2 /Ni = N R (N 2 ) 3 /N R (Ni) 2 arc universal in the sense 
that they are insensitive to structural details of polymer 
models. From such a consideration, these quantities are 
estimated as N R (Ni) ~ 6 and N R (N 2 ) ~ 15, which then 
yield g* ~ 6Ni and N 2 ~ 10 2 iVi. 

We thus find two well-defined regimes (i) the ideal 
statistic regime R = HN/g) 1 / 2 = a^^N 1 / 2 for short 
rings (N < N\) and (ii) the compact statistic regime 
R = TW.9*) 1/3 = a^/V 1 / 6 ^ 1 / 3 for long rings 
(N > N 2 ), which are connected through a rather wide 
crossover regime. By comparing the above two expres- 
sions for the ring size in respective regimes, one finds 
that they match at N = g*. However, these regimes (i) 
and (ii) are restricted in the range N < Ni — g*<f>% 
and N > N 2 = g* N R (N 2 ), a straightforward conse- 
quence of which is an effective exponent v e ff (if de- 
fined) in this crossover regime should lie in the range 
u e ff € (1/2,1/3). The closer analysis shows that the 
ring size obeys R = ^{N/g^) 1 / 3 = a^ -1 / 8 ^ 6 N 1 / 3 in 
this regime, with the scale-dependent mesh size £,^(N) = 



acf) 1/>s [g^ (N)] 1 / 2 . Such a feature should be manifested in 
the intra-molecular correlation (form factor) of individ- 
ual rings in melts. The evolution of ^(N) =C<t> R q) with 
N is controlled by the balance between inter-molecular 
and intra-molecular topological constraints. Its growth 
rate is indeed slow, limited by the logarithmic diver- 
gence in Fi n ter {4>r) with <j) R — > 1, and well fitted by 
the effective power law ^ ~ TV 1 /' 5 g^ ~ TV 2 /' 5 , yielding 
u e ff ~ 1/3+ (2/5)(l/6) = 2/5. Our theory thus claims 
that the exponent proposed in C-D conjecture is a fictivc 
one. Nevertheless, the presence of the wide crossover 
range put forward the practical usage of this fictive ex- 
ponent, based upon which various effective power laws 
on statics and dynamics in the crossover should be con- 
structed. 

We have discussed the underlying physical picture be- 
hind the excluded- volume analogy in Sec. IIII CI Here it 
should be emphasized that compared to the usual hard 
sphere systems, the unique property of ring polymer 
melts show up owing to their great flexibility and in- 
ternal degrees of freedom, which is incorporated in the 
present model through the Y-factor. This softness fac- 
tor is a consequence of the presence of intrinsic topo- 
logical length scale Ni, allowing partial overlapping (in- 
terpenetration) of different rings, and provides a leeway 
to form the topologically active network of meshes, i.e., 
</> R (N) = N R (N)Y < 1 at JVi < N < N 2 (eq. ©). A 
fictive exponent v e ff — 2/5, which is slightly larger than 
that for the compact statistics, results from the progres- 
sive interpenetration of surrounding rings. From the ob- 
servation that such a leeway is gradually lost with the 
increase in N R {N) towards N R (N 2 ), we find that the 
presence of this broad crossover regime is inextricably 
linked to the presence of the intrinsic length scale N±. 
Beyond the onset of the compact statistics N > N 2 , each 
ring resembles a hard ellipsoid in a statistical sense with 
the coordination number N R (N 2 ) equal to that in hard 
ellipsoid systems, but there should be, of course, sub- 
stantial size and shape fluctuations, and the rings can be 
rather intensively deformed by strong external forcing. 

The stiffness effect apparent in the scmiflexiblc rings is 
shown to be absorbed in N\. Here, it was demonstrated 
that the stiffness dependence of Ni (cq. (|35|) ) closely re- 
sembles that of N e in the corresponding linear chain melts 
(Fig|6|). This result together with our geometrical inter- 
pretation of iVi may provide an intriguing route towards 
a better understanding of the entanglement concept. 

In closing, we would like to emphasize that the melt of 
ring polymers is not a mere theoretical construction. In 
addition to past experimental attempts [ToL [l8l l37l [38| , 
rece nt prog resses in synthesis and characterization meth- 
ods [39144 1| are worth noting, which may open up the 
quest for unique material properties of topological ori- 
gin Q. Its close connection to the chromosome territo- 
ries, i.e., a higher-order spatial organization of chromo- 
somes inside the cell nucleus, provides a further motiva- 
tion for the investigation [3 a . |42j . There are also other 
related issues such as topological effects in phase scpara- 
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tion [43[ , confined ring polymer melts 44 1 and the prop- 
erty of brush composed of ring polymers (4f| . We hope 
that the present framework and the emergent geometri- 
cal picture will improve our understanding of polymeric 
systems under topological constraints. The elucidation 
of the dynamics and rheology [Io|, [HI Hlj of ring polymer 
melts from the present viewpoint might be an interesting 
challenge. 



Acknowledgments 



T.S thanks A. Takano for discussion on the experi- 
mental situation of ring polymer systems including sam- 
ple preparation, characterization methods and physical 
measurements. He is also grateful to J. Juul for his kind 
help in proofreading. 



[1] D.M. Raymer and D.E. Smith, Proc. Natl. Acad. Sci. 
USA 104, 16432 (2007). 

[2] P.-G. de Gennes, Scaling Concepts in Polymer Physics 
(Cornell University Press, Ithaca, 1979). 

[3] M. Doi and S.F. Edwards, The Theory of Polymer Dy- 
namics (Oxford University Press, Oxford, 1986). 

[4] M. Rubinstein and R.H. Colby, Polymer Physics (Oxford 
University Press, Oxford, 2003). 

[5] S.F. Edwards, Proc. Phys. Soc. 92, 9 (1967). 

[6] P.-G. de Gennes, J. Chem. Phys. 55, 572 (1971). 

[7] M.E. Cates and J.M. Deutsch, J. Physique 47, 2121 
(1986). 

[8] T. McLeish, Science 297, 2005 (2002). 

[9] A. Grosberg and A. Khokhlov, Statistical Physics of 

Macromolecules (AIP, NY, 1994). 
[10] M. Kapnistos et. al. Nature Mater. 7, 997 (2008). 
[11] T. Vettorel, A.Yu. Grosberg and K. Kremer Phys. Biol. 

6, 025013 (2009). 
[12] J.D. Halversonet. al. J. Chem. Phys. 134, 204904 (2011); 

ibid 134, 204905 (2011). 
[13] T. Sakaue, Phys. Rev. Lett. 106, 167802 (2011). 
[14] R. Everaers, et. al., Science 303, 823 (2004). 
[15] M.Miiller, J.P. Wittmer and M.E. Cates, Phys. Rev. E 

61, 4078 (2000). 
[16] M.Miiller, J.P. Wittmer and J.-L. Barrat, Europhys. 

Lett. 52, 406 (2000). 
[17] J. Suzuki, A. Takano, T. Deguchi and Y. Matsushita, J. 

Chem. Phys., 131, 144902 (2009). 
[18] A. Takano, Polym. Prepr. Jpn. 56, 2424 (2007). 
[19] M.D. Frank-Kamenetskii, A.V. Lukashin and A. V. Vol- 

ogodskii, Nature 258, 398 (1975). 
[20] J.des Cloizeaux, J. Phys. Lett. 42, L433 (1981). 
[21] T. Sakaue and E. Raphael, Macromolecules 39, 2621 

(2006). 

[22] E. Flikkema and G. ten Brinke, J. Chem. Phys. 113, 
11393 (2000). 

[23] P. Chaikin, A. Donev, W. Man, F. Stillinger and S. 
Torquato, Ind. Eng. Chem. Res. 45, 6960 (2006). 



[24] A.Yu. Grosberg, S.K. Nechaev and E.I. Shakhnovich, J. 

Phys. France 49, 2095 (1988). 
[25] T. Pakula and S. Geyler, Macromolecules 21, 1665 

(1988). 

[26] M.Miiller, J.P. Wittmer and M.E. Cates, Phys. Rev. E 

53, 5063 (1996). 
[27] S. Brown, G. Szamel, J. Chem. Phys. 109, 6184 (1998). 
[28] M. Daoud and J.F. Joanny, J. Physique 42, 1359 (1981). 
[29] A.R. Khokhlov and S.K. Nechaev, Phys. Lett. A 112, 

156 (1985). 

[30] S.P Obukhov, M. Rubinstein and T. Duke, Phys. Rev. 

Lett. 73, 1263 (1994). 
[31] ST. Milner and J.D. Newhall, Phys. Rev. Lett. 105, 

208302 (2010). 

[32] B.H. Zimm and W.H. Stockmayer, J. Chem. Phys. 17, 
1301 (1949). 

[33] T.C. Lubensky and J. Isaacson, Phys. Rev. A 20, 2130 
(1979). 

[34] D.W. Schaefer, J.F. Joanny and P. Pincus, Macro- 
molecules 13, 1280 (1980). 

[35] T. Sakaue, Macromolecules 40, 5206 (2007). 

[36] A.Yu. Grosberg, Y. Rabin, S. Havlin and A. Neer, Euro- 
phys. Lett. 23, 373 (1993). 

[37] J. Roovers, Macromolecules 21, 1517 (1998). 

[38] V. Arrighi, et. al., Macromolecules 37, 8057 (2004). 

[39] C.W. Bielawski, D. Benitez, R.H. Grubbs, Seience 297, 
2041 (2002). 

[40] A. Takano et. al., Macromolecules 40, 679 (2007). 

[41] Y. Ohta, Y. Kushida, Y. Matsushita and A. Takano, 
Polymer 50, 1297 (2009). 

[42] E. Lieberman-Aiden, et. al., Science 326, 289 (2009). 

[43] A.R. Khokhlov and S.K. Nechaev, J. Phys. II France 6, 
1547 (1996). 

[44] Private communication with Y.-L. Chen. 

[45] D. Reith, A. Milchev, P. Virnau and K. Binder, Euro- 
phys. Lett. 95, 28003 (2011). 



